Accepted for publication in The Astrophysical Journal 
Preprint typeset using IAT^X style cmulatcapj v. 10/09/06 



VARIABLE SODIUM ABSORPTION IN A LOW-EXTINCTION TYPE IA SUPERNOVA 1 - 2 

Joshua D. Simon 3 , Avishay Gal- Yam 4 , Orly Gnat 5,6 , Robert M. Quimby 7 , Mohan Ganeshalingam 8 , Jeffrey M. 
Silverman 8 , Stephane Blondin 9 , Weidong Li 8 , Alexei V. Filippenko 8 , J. Craig Wheeler 10 , Robert P. Kirshner 11 , 
Ferdinando Patat 9 , Peter Nugent 12 , Ryan J. Foley 1113 , Steven S. Vogt 14 , R. Paul Butler 15 , Kathryn M. G. 
Peek 8 , Erik Rosolowsky 16 , Gregory J. Herczeg 17 , Daniel N. Sauer 18 , and Paolo A. Mazzali 19 ' 20 ' 21 

Accepted for publication in The Astrophysical Journal 



o 
o 



X 

6 
to 



> 
m 
oo 
o 

o 

OS 

o 



X 



ABSTRACT 

Recent observations have revealed that some Type la supernovae exhibit narrow, time-variable 
Na I D absorption features. The origin of the absorbing material is controversial, but it may suggest 
the presence of circumstellar gas in the progenitor system prior to the explosion, with significant 
implications for the nature of the supernova progenitors. We present the third detection of such 
variable absorption, based on six epochs of high-resolution spectroscopy of the Type la supernova 
SN 20071e from the Keck I telescope and the Hobby-Eberly Telescope. The data span a time frame of 
approximately three months, from 5 days before maximum light to 90 days after maximum. We find 
that one component of the Na I D absorption lines strengthened significantly with time, indicating 
a total column density increase of ~ 2.5 x 10 12 cm -2 . The data limit the typical timescale for the 
variability to be more than 2 days but less than 10 days. The changes appear to be most promi- 
nent after maximum light rather than at earlier times when the ultraviolet flux from the supernova 
peaks. As with SN 2006X, we detect no change in the Ca II H and K absorption lines over the same 
time period, rendering line-of-sight effects improbable and suggesting a circumstellar origin for the 
absorbing material. Unlike the previous two supernovae exhibiting variable absorption, SN 20071e is 
not highly reddened (Eb-v = 0.27 mag), also pointing toward circumstellar rather than interstellar 
absorption. Photoionization calculations show that the data are consistent with a dense (10 7 cm~ 3 ) 
cloud or clouds of gas located ~ 0.1 pc (3 x 10 17 cm) from the explosion. These results broadly sup- 
port the single-degenerate scenario previously proposed to explain the variable absorption, with mass 
loss from a nondegenerate companion star responsible for providing the circumstellar gas. We also 
present possible evidence for narrow Ha emission associated with the supernova, which will require 
deep imaging and spectroscopy at late times to confirm. 

Subject headings: circumstellar matter — supernovae: general — supernovae: individual (SN 1999cl, 
SN 2006X, SN 20071e) 
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1. INTRODUCTION 

The origin of Type la supernovae (SNe la) has impor- 
tant implications for understanding the physics of these 
dramatic explosions, the evolution of binary stars, and 
the expansion of the universe. Recent interest in this 
last issue is particularly high because SNe la are the 
most effective distance indicators known on cosmological 
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scales. These objects are thus the focus of many current 
and future dark-energy experiments (e.g.. lAstier et al 



2006 [lRiess et al.ll2007tlWood-Vasev et al.l2007l:lAldering 
20051) . Investigating systematics that could affect their 
luminosity evolution in any way is a prerequisite for pre- 
cision cosmological studies using SNe la as calibratable 
standard candles. 

Perhaps the most serious gap in our understanding 
of SNe la is our lack of knowledge about their pro- 
genitor systems. It is widely thought that SNe l a 
occur in close binary systems ()Whelan fc Ibenl Il973f ). 
in which a carbon-oxygen white dwar f accretes mat- 
ter fro m another star until it nears the [Chandra sekharl 
<fl93lh mass and its core ignites a thermonuclear explo- 
sion that propagates outward and destroys the white 
dwarf. The identity of the companion star, however, 
is unknown; possibilities range from a second white 
dwarf (a double-de generate system; llben fc Tutukovl 
119841 : Wcbbink ll984[) to a main-sequence star, subgiant, 
or evolved red giant (all single-degenerate systems). At 
best, these proposed progenitor systems would have an 
absolute magnitude of My w —3 (for a red giant com- 
panion star near the tip of the red giant branch), leav- 
ing them undetectable by current telescopes for all but 
the closest galaxies. Si nce the last known Local Group 
SN la occurred in 1885 (|de Vaucouleurs fc Corwinlll985l ; 
iFesen. Saken. fc Hamilton 1989) and there have been no 
SNe la observed in t he Milky Way in over 400 years 
(iRuiz-Lapuentd 12004 [Badenes et all 120061 : iKrause et all 
2008), direct empirical constraints on the nature of the 
progenitors are severely lacking. 

Recent observations, however, have yielded progress 
on this issue using several different approaches. First, 
statistical studies of large numbers of SN host galaxies 
have demonstrated a strong relationship between Hubble 
type and the rate of Type la supernovae, with late-type 
galaxies hosting as many as 20 time s more SNe la per 
unit mass as early-type galaxies (e.g.. | Oemler fc Tinslevl 
119791: iMannucci etaU 120051 : ISullivan et al.ll2006h . This 
result suggests that many SNe la come from a rela- 
tively young stellar population, although of course the 
SNe la in elliptical galaxies must have their origin in 
an older population. These observations, combined with 
the SN la rate as a function of redshift, can be used to 
derive the delay-time distribution (the amount of time 
between the formation of a stellar population and the 
occurrence of SNe la). The measured delay-time dis- 
tribution includes both a prompt (~ 10 8 yr) compo- 
nent that is closely associated with recent star formation 
and a much more extended com ponent with a character- 
istic timescale of several Gvr (|Scannapieco fc Bildstenl 
120051: IMannucci Delia Valle fc Panagial l2006t although 
see lTotani et al.ll2008l for a conflicting view). These two 
distinct components of the SN la population would seem 
to naturally correspond to different mechanisms for pro- 
ducing Type la supernovae, and perhaps different pro- 
genitor systems as well. Supporting this idea, a number 
of studies have demonstrated that SNe la in early-type 
and late- type galaxies have different average luminosi- 
ties (e.g..lHamuv et al.l fl996l: lHowelll[200l ILi et al.ll200ll : 
iHicken et al.ll2009h . 

Detailed studies of indiv idual nearby SNe la have 
also re vealed crucial clues. iQuimbv. Hoflich. fc Wheeler! 
(|2007t ) found unusual behavior of the velocity of the 



Si II A6355 line in the otherwise normal SN 2005hj and 
argued that this object might be t he prototype of a sep- 
arate class of SN la explosions. iRuiz-Lapuente et al.l 
(|2004f ) obtained spectra of stars near the center of 
the remnant of Tycho's supernova (SN 1572) and 
claimed to identify the companion star to the SN pro- 
genitor, a G-type subgiant, althoug h this interpreta- 
tion has proved controversial (e.g.. Ilhara et al.l 120071 : 
iGonzalez- Hernandez et al.ll2009l : [Kerzendorf et al.ll2009h . 

Finally, another recent advance has been made by 
iPatat etail (|2007aD . who detected circumstellar mate- 
rial (CSM) in a SN la via optical absorption lines for 
the first time. Previous searches for emission from 
CSM in SNe la using Ha, radio, and X-ra y observations 
have generally yielded only upper limits (iMattila et al 
2005[ iLeonaTdl I2007t iPanagia et al] 120061 : llmmler et al 
2006 : iHughes et al.l 120071 ). with the exception of a 



few unusual objects whose classification as SNe la has 
been disputed (e.g., ILi et alJ l200ll 120031: lHamuv et al.1 
[ 20031: lAldering et a l. 2006: iPrieto et all I2007L but see 
iBenetti et al.1 [20061 and iTrundle et al.l 120081 for alterna- 
tive interpretations). By obtaining high-resolutio n spec- 
tra of SN 2006X at multiple epochs, IPatat et"aTI showed 
that the strength of at least four distinct components 
of the Na I D absorption lines varied with time. The 
temporal coverage of the data was rather sparse, but the 
variability occurred over the time span between 2 days 
before maximum light and 61 days after. Similar rapid 
changes in metal absorption-line profiles with time are 
commonly seen in novae and prov ide clues to the mass - 
transfer process in those systems (jWilliams et aLl feOOS). 
Variable Na D absorption is also observed in the spec- 
tra of Milky Way stars (generally on timescales of years 
to decades, although observations on shorter timescales 
are lacking) and is usually attributed to small interstel- 
lar clouds moving perpend icular to the line of sight (e.g., 
IWelty fc Fitzpatrickl l2001h . The lack of corresponding 
variability in the Ca II H & K absorption features in 
SN 2006X, however, renders this interpretation viable 
only with an appeal to peculiar chemical abundances 
and/or an unusual geometry in the absorbing interstellar 
clouds. 

Accordingly, IPatat et"aT1 (|2007al ) concluded that the 
variable absorption features originate in circumstellar 
clouds in the progenitor system that were ionized by the 
radiation from the supernova and recombined over the 
following several weeks, with some of the clouds then be- 
ing collisionally reionized by the SN ejecta at later times. 
This circumstellar gas could have originated either in the 
stellar wind from the progenitor's companion star or in 
successive nova eruptions. Because Na I has a much lower 
ionization potential than Ca II, the Na D line profiles 
can change without an accompanying effect in the Ca H 
& K lines if t he ionizing r adiation has an appropriate 
spectrum. The IPatat et al.l interpretation has been chal- 
lenged by IChugail (|2008f ). who showed that for typical 
red giant w ind den s ities N a D absorption should not be 
detectable. [Chugai (2008]) suggested instead that the ab- 
sorbing material must be located farther away from the 
supernova, possibly not associated with the progenitor 
system at all. Photometry and late-time low-resolution 
spectroscopy of SN 2006X, which reveal at least one 
light echo from the supernova, als o provide support f or 
a dusty circumstellar environment (jWang et al.ll2008aU bl: 
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ICrotts fe Yourdonll2008D . Taken as a whole, these re- 
sults appear to indicate a single-degenerate progenitor 
for SN 2006X, perhaps with a red giant companion. 

But is SN 2006X merely a unique, pathological object? 
Or is it broadly represen tative of a sign if icant fraction of 
(or perhaps all) SNe la? iBlondin et a l. (2009) have now 
identified a second supernova, SN 1999cl, displaying vari- 
able Na D absorption that is visible even at low spectral 
resolution, but the lack of high-resolution data or cov- 
erage of the Ca H & K lines m akes the interp retation 
of this object more difficult. The IBlondin et all analysis 
demonstrates that significant changes in the Na D ab- 
sorption profile are relatively uncommon, with only two 
detections out of 31 SNe examined. Notably, SN 1999cl 
and SN 2006 X are the two most h eavily reddened super- 
novae in the IBlondin et alj (|2009f) sample, creating the 
appearance of a connection between large (presumably 
interstellar) reddening and variable absorption. 

Only two other SNe la (SN 2000cx and SN 2007af) 
have multi-epoch high-resolution spectra available in the 
literature, and neither of those objects e xhibits any 
changes in the Na D line profiles with time (|Patat et al.l 
2007b!; ISimon et al.ll2007T ). In this paper we present a 
similar data set for the bright SN la 20071e, including six 
high-resolution spectra obtained between —5 and +90 
days rel ative to maximum light. We use these data to 
test the lPatat et al.l (j2007al ) scenario, searching for vari- 
ability in the Na D absorption features. We describe all 
of our observations, both photometric and spectroscopic, 
in § [2j In § [3l we compare the light curve and spectral 
evolution of SN 20071e with those of other SNe la and 
analyze the high-resolution spectra. We discuss the im- 
plications of the variable Na D absorption features in § [4] 
and present our conclusions in § O 

2. OBSERVATIONS AND DATA REDUCTION 

SN 20071e was discov e red by L. Monard 
(|Monard. Yamaoka. fc Itagakil I2007D on 2007 Octo- 
ber 13.79 (UT dates are used throughout this paper). 
Optical spectra obtained two nights later showed that 
the object was a Type la supernova at least one week 
before maximum light, and featurin g high-velocity 
ejecta expanding at ^16,000 km s _1 (Filippcnk o et alJ 
120071) . The host galaxy of the supernova is NGC 7721, 
an Sc galaxy with a heliocentric r ecession velocity of 
2015 km s" 1 (jKoribalski et al.ll2004D . 

2.1. High- Resolution Spectroscopy 

Our high-resolution observing campaign for SN 20071e 
began on 2007 October 20 w ith the High Re solution 
Echelle Spectrometer (HIRES; iVogt et al.lll994f) on the 
Keck I 10 m telescope. Over the following three months, 
we obtained a total of five HIRES spectra of the su- 
pernova. The seeing was often poor (> 1") during the 
observations, and the overall conditions varied signifi- 
cantly over the course of the many observing runs. The 
data were obtained with a range of different spectro- 
graph setups, including spectra with both the blue and 
red cross-dispersers, and achieved signal-to-noise ratios 
(S/N) ranging between 11 and 113 per pixel. Most of 
the spectra used a 7.0" x 0.86" slit, yielding a spectral 
resolution of R ss 52, 000, but one spectrum was obtained 
with a wider (7.0" x 1.15") slit and R « 42,000. All of 
the HIRES spectra cover the Ca H & K, Na D, and Ha 



lines. We reduced the Keck data using the IDL data re- 
duction package for HIRES (version 2.0) developed by 
J. X. Prochaska and collaborators (Bernstein et al., in 
prep.). 22 

We also obtained one spectrum of SN 20071e w ith the 
High-Resolution Spectrograph (HRS: lTulllll998l ) on the 
Hobby-Eberly Telescope (HET) on 2007 November 4. 
The spectrograph was in its R = 60, 000 mode, with a 
2"-diameter fiber and the 316 line/mm grating centered 
at 6948 A, providing nearly complete wavelength cover- 
age from 5100 to 8850 A. We obtained 3 spectra totaling 
3000 s of exposure time and reached a combined S/N 
of 16 per pixel. The HRS data were reduced in IRAF 23 
with the echelle package using standard procedures. A 
comprehensive summary of all the high-resolution data 
is presented in Table [T] 

2.2. Imaging and Low-Resolution Spectroscopy 

SN 20071e was the target of extensive photomet- 
ric follow-up observations with the 0.7 6 m Katzman 
Automatic Imaging T elescope (KAIT; iLi et al.l 120001: 
iFilippenko et al.l l200lD . continuing for approximately 3 
months until the supernova went into conjunction with 
the Sun. We obtained post-explosion images in 2008 Au- 
gust after the SN had faded in order to subtract t he host- 
galax y light. We used the daophot package ([Stetson! 
1987) in IRAF to perform point-spread function pho- 
tometry of SN 20071e relative to various field stars in the 
KAIT images, which were calibrated on 5 photometric 
nights with KAIT and the Nickel 1 m telescope at Lick 
Observatory. 

We also obtained low-resolution spectra of SN 20071e 
with the Low R esolution Imaging Spectrometer (LRIS; 
lOke et aT1ll995l ) on the Keck I telescope on 2007 October 
15, October 16, November 11, November 12, and Decem- 
ber 1 2, and with the Kast spectrograph (Mill er fc Stond 
1993) on the 3 m Shane telescope at Lick Observa- 
tory on 2007 November 2, November 18, and December 
1. These data were reduced in IRAF and IDL follow- 
ing normal procedures (for details , see lFolev et al.ll2003l : 
iMatheson etalll2000t lHorndll986l ). 

3. DATA ANALYSIS AND RESULTS 

3.1. Light Curves and Low- Resolution Spectra 

We display BVRI light curves of SN 20071e in Figure 
[TJ We fitted the photometric data with the latest version 
of the multicolor light-curve shape method (MLCS2k2; 
Uha. Riess. fc Kirshnej[2007l) to determine the parame- 
ters of the supernova. We find that the time of -B-band 
maximum was 2007 October 25.65 (JD = 2,454,399.15), 
with an uncertainty of 0.09 day. The derived line-of- 
sight extinction to the SN is Ay — 0.71 ± 0.06 mag, 
with an extinction law of Ry = 2.56 ± 0.22 (the Milky 
Way foreground reddeni ng is 0.033 mag, correspond- 
ing t o Ay = 0.11 mag; ISchlegel. Finkbeiner. fc Davis! 
Il998h . The distance modulus to SN 20071e is m - M = 
(32.35 - 51og[H /(70 km s" 1 Mpc" 1 )]) ± 0.06 mag, giv- 

22 Documentation and code for this reduction package are avail- 
able at http://www.ucolick.org/~xavier/HIRedux/index.html 

23 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 1 

High-Resolution Spectroscopy Observing Log 



Telescope Instrument 


UT Date 


Supernova Epoch a Exposure Time S/N b Spectral Resolution 


Wavelength Range 






[b] (A/AA) 


[A] 



Keck I 


HIRESb 


2007 October 20.36 


day —5 


7200 


113 


52,000 


3840- 


-6733 


Keck I 


HIRESr 


2007 October 25.34 


day 


900 


71 


41,700 


3460- 


-7989 


HET 


HRS 


2007 November 4.15 


day +10 


3000 


16 


60,000 


5100- 


-8850 


Keck I 


HIRESb 


2007 November 6.39 


day +12 


4500 


76 


47,600 


3845- 


-6649 


Keck I 


HIRESr 


2008 January 17.23 


day +84 


3300 


16 


54,000 


3870- 


-8365 


Keck I 


HIRESr 


2008 January 23.21 


day +90 


1200 


11 


54,000 


3900- 


-8360 



a Relative to maximum light. b Signal-to -noise ratio measurements are per pixel at the wavelength of the redshifted Na D lines (5934.5-5936.5 A). c The 
spectral resolution was determined by measuring the FWHM of emission lines in the comparison-lamp spectra and telluric absorption lines in the SN 
and telluric standard star spectra. 




Fig. 1— BVRI light curves (filled symbols) of SN 20071c, 
from KAIT and Lick Nickel 1 m data. Photometric uncertain- 
ties are indicated by the plotted error bars, which in most cases 
are smaller than the displayed data points. The solid lines repre- 
sent the MLCS2k2 fits to the data and the shaded bands represent 
the uncertainties on the fits. 

ing the SN an absolute magnitude of My = (—19.34 + 
5 log Ho/70)±0.09. The luminosity /light-curve-shape pa- 
rameter is A = -0.14 ±0.02, and the MLCS2k2 \ 2 value 
of 163.0 for 106 degrees of freedom indicates an accept- 
able fit (much of the x 2 comes from the /-band data, 
which significantly differ from the fit). 

In Figure [2] we show the spectrum of SN 20071e at 
an epoch of 10.3 days before maximum li ght. Analysis 
with t he Superfit spectral fitting code of iHowell et al.1 
(2005) indicates that SN 20071e is most similar to spec- 
tra of the highly reddened type la SN 2 002bo at 11 
days before maxim um from iBenetti et ail (|2004[ ). while 
the SNID package (jBlondin fc Tonrvl 120071 ) finds a best 
fit to early spectra of the "g olden standard" type la 
SN 2005cf (|Wang et al.l f2009h . Spectra of SN 2002bo 
and SN 2005cf were included in the databases used by 
both fitting packages; the differences between the two lie 
in their f itting methods and trea tment of color informa- 
tion (see iBlondin fe Tonrvl l2007f ). but since SN 2002bo 
and SN 2005cf were similar events in many respects 
(e.g., line velocities) we do not regard these fit re- 
sults as a signi ficant disagreemen t. Like SN 2002bo 
(and SN 2006X: IWan"g et al.l[2008ah . SN 20071e exhibits 
a high-velocity component to the Ca II near-infrared 
triplet absorption lines, although such features are very 



Fig. 2. — Low-resolution spectrum of SN 20071e 10.3 days be- 
fore maximum light (lower cur ve), compared to scaled spectra of 
SN 2005cf (mid dle, data fro m Wang et al. 2009]) and SN 2002bo 
(top, data from IBenetti et al.H2004l ) at similar epochs. The spec- 
tral features match closely between the three supernovae. The 
SN 2005c f and SN 2002bo data have been dere ddened by Ay = 
0.58 mag l lWang et al.H2009fl and A v = 1.0 mag (Elias-Ro sa et al.l 
2008), respectively, and then re-reddened by Ay = 0.71 mag to 
match the extinction of SN 20071e. 



common in SNe la that are obs erved at sufficiently 
early times (|Mazzali et al.l l2005bO . Both SN 2006X 
and SN 2002bo were highly p olarized (|Wang et al.ll2006l . 
I2007t IWang fc Whe eler 2008), suggesting that polariza- 
tion measurements of SN 20071e might be very interesting 
as well. 

In order to classify SN 20071e among the various known 
subgroups of Type la supernovae, we analyzed th e spec - 
tra according to the pres c riptio ns of lBenetti et al.l (2005) 
and iBranch et ail (|2006l l2009t) . The widths of the ab- 
sorption features near 5750 A a nd 6100 A co rrespond 
to the broad line (BL) group of IBranch et al.l . but the 
equivalent widths are not too far from those of the core 
normal group. We measured a velocity gradient in the 
Si II A6355 line of 83±3 km s" 1 day" 1 , placing SN 20071c 
m arginally with i n the high velocity gradient (HVG) class 
of lBenetti et alj (|2005f ). although this gradient is not very 
much larger than is seen in some low velocity gradient 
SNe. SN 2006X and SN 2002bo both exhibited signifi- 
cantly larger gradients. 

3.2. Removal of Telluric Absorption Features 
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The host-galaxy Na D absorption lines of SN 20071e 
are located at observed wavelengths between 5930 and 
5940 A. This region of the spectrum unfortunately con- 
tains a number o f telluric absorption lines from wa - 
ter molecules fe.g.. lMoore. Minnaert. &: Houtgastlll966f h 
Before searching for variations in the intrinsic super- 
nova line profiles, we therefore must remove the telluric 
featu res, which will change with atmospheric con ditions 
(e.g. JWade fc HorndH^a iMatheson et al.ll2001l ). 

On every night except that of the second-epoch obser- 
vations (day 0, 2007 October 25) we obtained at least 
one spectrum of a hot, rapidly rotating star to serve 
as a telluric standard. These spectra were reduced and 
normalized to a flat continuum level in the same man- 
ner as the supernova data. Because the telluric stan- 
dards do not contain any atmospheric features in the 
region of interest, any absorption seen must necessar- 
ily be of telluric origin. From the telluric standard ob- 
servations we constructed telluric absorption line model 
spectra using a list of known telluric lines. 24 Pixels in 
the model spectra that were at the wavelength of telluric 
lines were set equal to the value of those pixels in the tel- 
luric standard star spectra, and all other pixels were set 
to unity. We then selected the telluric model obtained 
at the most similar airmass and time for each supernova 
spectrum. 25 The telluric models were further adjusted 
to match the airmass of the supernova observations with 
the following scaling: model a djusted = (model or i gi „ a ;) b , 
where b = airmasssN/dirmass s tandard- Finally, we di- 
vided the reduced supernova spectra by their respective 
telluric absorption models to produce clean SN spectra. 
The strength of the absorption lines that overlap with 
the wavelength of the host-galaxy Na D absorption in 
the telluric models ranges from a few percent to 10%, 
not nearly large enough to account for the absorption 
profile changes described in §3.31 

3.3. Sodium D Absorption Lines at High Resolution 

We display the high-resolution spectra of SN 20071e 
around the host-galaxy Na D lines in Figure [31 At least 
one component of the absorption (near 5931.5 A and 
5937.5 A in the D 2 and Di lines, respectively) strength- 
ens over the course of the observations, increasing from a 
total depth of ~ 0.6 in the first spectrum to nearly satu- 
rated three months later. 26 The full absorption profile is 

24 We use a model rather than simply the observed standard star 
spectrum in order to remove any effects from imperfectly fitting 
the continuum of the standard star and because the S/N of the 
standard spectra is in some cases not too much higher than the 
S/N of the supernova spectrum. 

25 For the second-epoch observations where we lack a telluric 
standard from the same night, we use a telluric standard from the 
night of 2007 October 20, which had the most similar atmospheric 
water vapor content to 2007 October 25 according to the Mauna 
Kea Weather Center forecasts (available on the internet) for the 
nights on which wc observed. 

26 In Figure [3] it appears that the final (day +90) spectrum 
shows additional changes, with all of the Na absorption features 
weakening compared to the previous epochs. Given the complete 
lack of variations in all but one of these components through the 
day +84 spectrum obtained just six days earlier, the very poor ob- 
serving conditions under which the day +90 spectrum was acquired 
(bad seeing, high humidity, and bright sky because the supernova 
set before the end of twilight), and the low S/N of the spectrum, 
we suspect that these changes indicate a systematic problem with 
the data rather than true physical changes. We therefore exclude 
the day +90 spectrum from the remainder of our analysis. 



TABLE 2 
Na D Equivalent Widths 



Full Profile Red Side 



Epoch 


D 2 EW 
[mA] 


Di EW 

[mA] 


D 2 EW 
[mA] 


Di EW 
[mA] 


day —5 


894 + 3 


649 + 3 


653 + 2 


498 + 2 


day 


883 + 5 


661 ± 5 


655 + 3 


512 + 3 


day +10 


906 ± 20 


695 ± 20 


670 ± 12 


537 ± 12 


day +12 


949 + 4 


702 + 5 


711 + 3 


549 + 3 


day +84 


1006 ± 7 


766 ± 7 


759 + 4 


603 + 4 



complex, with ~ 7 distinct velocity components visible. 
Fitting the absorption with a set of Gaussians reveals 
that at least 2 additional components blended with the 
strongest absorption features are needed in order to re- 
move significant residuals, and even then the fits to the 
high-S/N spectra are not statistically satisfactory. Fits 
of multiple blended Gaussians do not produce unique re- 
sults, however, so we cannot study the physical condi- 
tions in the variable absorption component (s) accurately 
with this technique. 

Instead, we first integrate directly over the entire ab- 
sorption profile to determine the equivalent width (EW) 
of the Di and D 2 lines in each spectrum. The mea- 
sured equivalent widths are listed in columns (2) and 
(3) of Table[2j Over the course of the observations, both 
lines show an increase in the EW of slightly more than 
100 mA. The total EWs in the day —5 spectrum are 
894 ± 2 mA and 649 ± 2 mA for the D 2 and Di lines, 
respectively. By the final spectrum at day +84, the EWs 
are 1006 ± 5 mA and 766 ± 5 mA. For the four blue-most 
components of the absorption, we rule out any variation 
at high significance; the EWs of those components agree 
in all six epochs to a level of a few mA. To reduce the 
uncertainty in the estimate of the change in EW, we can 
therefore leave these components out of the integration 
and sum the absorption starting from the peak at 5931 A 
(5937 A for the Di line) instead of from the blue edge of 
the absorption at 5929.5 A (5935.5 A). With this smaller 
wavelength range we find a total increase in EW between 
day —5 and day +84 of 106 ± 5 mA and 105 ± 5 mA for 
the D 2 and Di lines (last two columns of Table [2|). 

It is worth noting that all of the Na D absorp- 
tion features, including the variable component, are 
blucshiftcd relative to the strongest absorption compo- 
nent at 5931.8 A (2130 km s _1 ), which we presume cor- 
responds to the local interstellar medium (ISM) velocity 
(see §313 and [3/7]). In SN 2006X and SN 1999cl, the 
variable absorpt ion was also bluesh i fted compared to th e 
ISM absorption (jPatat et al.ll2007al:lBlondin et al.ll2009f) . 
perhaps providing a clue as to the origin of the material 
responsible for the varying absorption. 

Having determined the EW of the variable absorption, 
we would now like to constrain the column density and 
line width of the absorbing gas. To isolate the changing 
component of the absorption profile, we examine the dif- 
ference between each spectrum and the first one. These 
difference spectra are displayed in Figure [4] At day 
(maximum light), we do not detect any obvious differ- 
ence signal from 5 days earlier, despite the possible small 
increase in EW reported for the Di line in Table [2j Be- 
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5928 5930 5932 5934 5936 5938 5940 5942 

Observed Wavelength (A) 



Fig. 3. — High-resolution observations of the Na D absorption lines in the spectrum of SN 20071e near the host-galaxy velocity. The top 
panel shows the initial HIRES spectrum (day —5) in black and the second-epoch spectrum from day later in red. In each subsequent 
panel the initial spectrum is plotted as a dotted black curve. The second panel compares the day (red) and day +10 (blue) spectra to 
each other, and then each following panel does the same for the next neighboring pair of observations. All of the spectra shown in each 
panel are smoothed to the lowest resolution spectrum displayed in that panel (which means that the effective resolution changes somewhat 
from panel to panel). At least one component of the Na absorption near 5931.5 A in the D2 line clearly strengthens with time. 



tween day —5 and day +10, however, a clear Gaussian 
line emerges in the difference spectrum (Fig. 01 second 
row from the top). As expected from the EW measure- 
ments above, the strength of this line increases through 
day +84, and we also find that the Doppler parameter 
of the line grows from 4.2 ± 0.7 km s _1 in the (day +10 
minus day —5) spectrum to 10.7 ± 0.2 km s _1 in the 
(day +84 minus day —5) spectrum (see Table [3]). The 
line in the (day +10 minus day —5) difference spectrum 
appears close to symmetric, but stronger hints of asym- 
metry in the later difference spectra suggest that some of 
this increase in line width is resulting from the addition 
of a second blended absorption component. 



3.4. Column Density of the Variable Component 

We see in Table [3] that the change in the EW of the 
D2 line is very similar to the change in the EW of 
the Di line at all epochs. At 1 a significance, 0.78 < 
AEW Di /AEWd 2 < 1.03 over th e full set of difference 
spectra. Following" ! Spitzerl (|1978l Equations 3-49 to 3- 
51), we can use a curve-of-growth analysis to determine 
the observed EW for any combination of column density 
and line width. Using the derived Doppler parameters 
for the difference spectra given in § 13.31 we calculate the 
change in EW that would be observed for a small in- 
crease in column density for all column densities between 
10 9 cm -2 and 10 14 5 cm -2 . For a given line width, these 
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TABLE 3 
Na D Difference Spectra Fits 



Epochs 


D 2 


EW change 


Di EW change 


D2 velocity D2 


Doppler parameter 






[inA] 


[mA] 


[km s _1 ] 


[km s" 1 ] 


day +10 — day 


-5 


50 ±4 


45 + 5 


2123.1 ± 0.2 


4.2 ± 0.7 


day +12 — day 


-5 


55+1 


54 + 2 


2122.9 ± 0.1 


5.3 ± 0.2 


day +84 - day 


-5 


85 + 2 


85 + 2 


2120.8 ± 0.1 


10.7 + 0.2 




5929 5930 5931 5932 5936 5937 5938 5939 

Observed Wavelength (A) Observed Wavelength (A) 
Fig. 4. — Differences in the Na D lines between later epochs 
and the first spectrum. The left panels show the D2 lines and the 
right panels the Di lines. The top curves in each panel are the 
observed spectra, and the black line at the bottom of each panel is 
the difference between the two spectra plotted in that panel. The 
orange curves are Gaussian fits to the difference spectra. There 
are no significant differences between the first two spectra (at and 
before maximum light), but the later spectra all show additional 
absorption at 5931.7 and 5937.7 A. 




10 11 12 13 

Log column density (cm -2 ) 

Fig. 5. — Ratio of the change in the EW of the Na I Di line to 
that of the D2 line for a small increase in the absorbing column 
as a function of column density. The dotted red curve represents 
the EW ratio for a Doppler parameter of b = 4.2 km s — 1 , appro- 
priate for the (day +10 minus day —5) spectrum, the dashed blue 
curve for a Doppler parameter of b = 5.3 km s _1 , appropriate for 
the (day +12 minus day —5) spectrum, and the dot-dashed gray 
curve for a Doppler parameter of b = 10.7 km s , appropriate 
for the (day +84 minus day —5) spectrum. The horizontal solid 
lines delineate the 1 <r constraints placed by the observed EW ra- 
tios, 0.78 < AEW Di /AEWd 2 < 1-03. The observed ratio of the 
EW changes can only be obtained over a narrow range of column 



densities from 10 



to 10 1 



sured line width larger than the true line width as noted 
in § 13.31 However, without a higher-resolution spectrum 
that separates the two components, we do not have con- 
clusive evidence in favor of this interpretation. 



calculations demonstrate that increasing the equivalent 
width of the D2 and Di lines by the same amount is only 
possible for a narrow range of column densities, as illus- 
trated in Figure O For b — 4.2 km s _1 , we find that 

12.0 < log7V Na 1 < 12.3, for b = 5.3 km s" 1 , we find that 

12.1 < logiV Na 1 < 12.4, and for b = 10.7 km s - \ we 
find that 12.4 < log N Na T < 12.7. We therefore conclude 
that the column density of the varying component of the 
absorption is A^a 1 ~ 2.5 x 10 12 cm~ 2 . Depending on the 
Doppler parameter, such a column density corresponds 
to a total Na D 2 EW of 200-300 mA for this component, 
consistent with the results shown in Table [2j 

It is evident from Figure [5] that the column densities 
indicated by the first two difference spectra (blue dashed 
and red dotted curves) are not consistent with the col- 
umn density preferred by the third difference spectrum 
(gray dot-dashed curve). We suspect that the reason for 
this inconsistency is that the third difference spectrum 
contains multiple blended components, making the mea- 



3.5. Calcium H&K Absorption Lines at High 
Resolution 

The Keck/HIRES spectra extend far enough to the 
blue that we also detect Ca II H & K absorption lines 
from both the Milky Way and the host galaxy, although 
the S/N of the day +90 spectrum in the blue is too low 
to be useful. We compare the host-galaxy absorption 
profile in Ca K and Na D 2 in Figure [H where a close 
correspondence between the absorption components in 
each species is visible. 

Unlike the Na lines, though, we detect no statisti- 
cally significant changes in the Ca absorption profile with 
time. Even nearly three months after maximum light, the 
profile shape (Figure [7| and total EW (Table [3| match 
those from before maximum light within the uncertain- 
ties. However, while changes in the Ca EW as large as 
those seen in the Na lines can be ruled out, the low S/N 
of the late-time spectra prevents us from being able to 
place strong constraints on smaller variations. 
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Heliocentric Velocity (km s -1 ) 

Fig. 6. — Comparison between the host-galaxy absorption profile 
in the Ca II K (upper spectrum, in red) and Na I D2 (lower spec- 
trum, in black) lines, both from five days before maximum light. 
The vertical dashed gray line is at the central velocity of the vari- 
able Na absorption component. Each Na absorption component, 
including the variable one, has a visible Ca counterpart (and vice 
versa), with the possible exception of some very weak blueshifted 
features between 1950 and 2000 km s — 1 . The relative strengths 
of the Na and Ca lines, however, vary significantly from one com- 
ponent to the next. The Ca lines also appear to have much more 
prominent wings than the Na absorption. The vertical ticks mark 
the velocities of individual absorption components inferred from 
Gaussian fits to the Na lines, and the vertical dotted magenta line 
indicates the velocity of the Hcv emission line discussed in § 13.71 



TABLE 4 
Ca H & K Equivalent Widths 





Full Profile 


Red Side 


Epoch 


K EW 
[mA] 


H EW 

[mA] 


K EW 
[mA] 


H EW 

[mA] 


day —5 
day 
day +12 
day +84 


653 ± 2 
677 ±9 
629 ± 6 
749 ± 32 


382 + 2 
397 + 9 
388 + 6 
413 ± 28 


301 ± 1 
308 + 4 
294 + 3 
343 ± 16 


174 + 1 
178 + 5 
186 + 3 
153 ± 15 



3.6. Other Absorption Features 

In addition to the Ca H & K and Na D lines, we also 
searched the spec tra for the addition al a bsorption fea- 
tures identified bv lPatat et alj (|2007aft and lCox fc Patatl 
dHH in SN 2006X. In the highest S/N (day -5) spec- 
trum, we detect the CH+ 3957.70 A and 4232.55 A lines 
at a velocity of 2130 km s -1 , exactly matching the veloc- 
ity of the strongest Na D absorption. Diffuse interstellar 
bands are also visible at 5780 A and 6283 A (rest wave- 
lengths). We do not detect Ca I at 4226.73 A, CH at 
4300.30 A, the CN vibrational band, or K I at 7699 A. 

3.7. Ha Emission 

All six of our high-resolution spectra cover the ex- 
pected wavelength of the redshifted Ha line. However, 
the fiber spectrum from the HET is not suitable for in- 
vestigating possible Ha emission from the SN because of 
the lack of local sky subtraction, so we only consider the 
five high-resolution Keck spectra in this section. 

We detect a narrow Ha emission line at the position 
of the supernova at all epochs of our observations. The 



HIRES spectra were obtained with a slit 7" long and ei- 
ther 0.86" or 1.15" wide, and faint Ha emission from the 
host galaxy is visible all the way along the slit. How- 
ever, careful sky subtraction reveals that there is also 
additional emission that is spatially unresolved and co- 
incident with the position of the supernova (see Figures 
[5] and . This line is also at a slightly different velocity 
than the extended host-galaxy emission. 

From the HIRES data, we measure the EW, veloc- 
ity, and line width of the Ha emission at the position 
of the supernova; the results are listed in Table [5] At 
v = 2139 km s _1 , the Ha emission is offset by 16 km s _1 
from the velocity of the variable sodium line (and has 
a higher velocity than any of the main sodium compo- 
nents), but is similar to the host-galaxy rotation velocity 
at the radius of the supernova measure d from a long- 
slit spectrum by lAfanasvev et al.l ([1992) . We then use 
the observed light curve and a comprehensive set of low- 
resolution spectra of SN 20071e from the CfA Supernova 
Archive to derive absolute flux calibrations for the spec- 
tra. In order to remove contaminating host-galaxy light 
(since the CfA spectra were obtained with a wide slit), 
the spectra are first scaled to the observed l/-band mag- 
nitudes of the supernova and then warped to match the 
B magnitudes as well. From these flux-calibrated spec- 
tra, we determine the SN continuum level at the posi- 
tion of Ha (6610 A), scaling the spectra up and down as 
necessary to account for the change in the supernova's 
R magnitude between the time of observation for each 
high-resolution spectrum and the nearest low-resolution 
spectrum (the time differences are less than two days for 
days —5 and +12, six days for day 0, and five and eleven 
days for days +84 and +90, respectively). We then con- 
vert the high-resolution EWs to fluxes in absolute units 
with the measured continuum levels. 

In principle, the Ha emission could originate cither 
from an H II region projected within 140 pc (correspond- 
ing to the «1" slit width) of the explosion or from the 
supernova itself. In the case of an H II region, one 
would expect other lines such as [O III] AA4959, 5007, 
[N II] AA6548, 6583, and [S II] AA6717, 6731 to be visible 
depending on the S/N. In the day +84 spectrum, which 
is the best for this purpose because of its relatively high 
S/N at red wavelengths and the faint magnitude of the 
SN at late times, we indeed detect H/3 and both lines 
of the [N II] and [S II] doublets. [O III], though, which 
is generally the brightest of the forbidden metal lines in 
H II regions (and comparable to or brighter than H/3), re- 
mains undetected at all epochs. [O III] fluxes more than 
an order of magnit ude weaker than Ha are uncommon in 
H II regions C e.g..lKennicutt. Bresolin. fc Garnet! 120031 : 
iMagrini et al.ll2007f ). 

On the other hand, if the Ha emission is related to 
the supernova, changes in the Ha flux with time would 
seem likely (although since such emission has never been 
detected, its exact characteristics are unknown). From 
the day —5 spectrum through the day +12 spectrum 
(spanning about 2.5 weeks), we detect no variability 
in the flux. In the day +84 spectrum more than two 
months later, however, the Ha flux is ~ 40% higher than 
at day +12. If the comparison is made only between 
days +12 and +84, the significance of this change is just 
2 a, but if we combine the measurements from days —5, 
0, and +12 together then the day +84 increase is signif- 
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5929 5930 5931 5932 5933 3959 3960 3961 3962 3963 

Observed Wavelength (A) Observed Wavelength (A) 

Fig. 7. — Comparison of the Na D2 (left panels) and Ca K (right panels) line profiles as a function of time. The top row includes spectra 
from day —5 and day 0, the middle row includes spectra from day and day +12, and the bottom row includes spectra from day +12 and 
day +84. The dotted black curve in the middle and lower panels represents the day —5 spectrum. Over the time period when significant 
changes are visible in the Na lines there are no detectable variations in the Ca lines. 



icant at the 4 a level. The low-S/N day +90 spectrum 
is consistent with both the early time and day +84 mea- 
surements. Taken as a whole, we consider these results to 
be a tentative indication of Ha emission from SN 20071c. 
More conclusive evidence of the nature of the Ha emis- 
sion awaits late-time imaging and spectroscopy to deter- 
mine whether or not there is an H II region coincident 
with the position of the supernova. 

If the Ha emission is associated with the SN, its lumi- 
nosity is related to the amount of gas present. The ob- 
served Ha emission has a luminosity of 1.6 x 10 35 ergs s _1 
(for the average Ha flux of 1.5 x 10~ 16 ergs cm" 2 s _1 ). 
If we follow the calculations of iLeonardl (|2007f ) , such 
a luminosity corresponds to ~ 0.025 M of solar- 
abundance material in the circumstcllar environment, 
in reasonable agreement with theoretical predictions 
dMarietta. Burrows & FrvxeH|[2000HMattila et al.ll2005l: 
iPakmor et al]l2008l ). However, those computations were 
specifically aimed at modeling Ha emission at very late 



times (day +380) and may not apply to the much ear- 
lier epochs at which we o bserved SN 20071e. Using the 
simple model presented bv lPatat et alj (|2007al ) as an al- 
ternative suggests a hydrogen mass of ~ 0.001 M Q for 
the observed Ha luminosity. The presence of hydrogen 
in the CSM has also been indirectly inferred from the 
broadening of absorption features in the earliest spec- 
tra of SNe la (jMazzali et alJ l2005a). In this scenario, 
~ 0.005 Mq of material with a solar composition is nec- 
essary to provide the electron density that is required to 
trap photons and ca use absorption in lines such as Ca H 
fc K or Si II A6355 (jMazzali et al.l l2005bb iTanaka et all 
120081) . 

4. DISCUSSION 

For SN 2006X. lPatat et all (|2007ah suggested a model 
in which the progenitor system consisted of a white dwarf 
and a red giant, with either the stellar wind from the red 
giant or repeated nova eruptions blowing off several pos- 
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Fig. 8. — Two-dimensional spectra of SN 20071e around the redshifted Ha line. These data are from day +84, where the S/N of the 
Ha detection is highest. The top panel shows a portion of the echelle order containing Ha in the coadded raw frames (the only processing 
was to remove cosmic rays). Diffuse Ha emission is visible across the slit, with a velocity that changes from one side to the other, but an 
additional distinct component is present right on top of the SN continuum. In the lower panel, the continuum emission and the average sky 
light across the slit have been subtracted. The sky fitting is not perfect, since the diffuse Ha was not symmetrically distributed around the 
supernova, but it is clear that there is excess Ha light at the position of the SN compared with the surrounding areas of the host galaxy. 
Fully reduced and extracted one-dimensional spectra of this emission are displayed in Figure [9] 



sibly asymmetric shells of material before the explosion. 
These circumstellar shells were proposed to be the lo- 
cation of the material responsible for the variable Na D 
absorption. Given timescales of the order of decades for 
the mass loss and the observed velocities of the variable 
lines, the absorbing material would be ~ 10 16 
from the star at the time of explosion 
envisions that the Na I atoms (with an ionization po- 
tential of 5.1 eV) are ionized by the UV radiation from 
the supernova, and then slowly recombine over the fol- 



cm away 
This scenario 



lowing weeks to produce strengthening Na D absorption 
lines. The absence of variations in the Ca H & K lines is 
then attributed to the much higher ionization potential 
of Ca II (11.9 eV), which prevents most of the Ca ions 
from being affected by the SN radiation field. 

However, IChugail (|2008[ ) has challenged this interpre- 
tation, arguing that the physical conditions expected in 
a red giant wind are not compatible with strong Na D 
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Fig. 9. — High-resolution spectra of SN 20071e around the red- 
shifted Ha and [N II] lines. From top to bottom, the spectra are 
from days —5, 0, +12, +84, and +90 (the day +10 spectrum is 
excluded because local sky subtraction cannot be done on a fiber 
spectrum). A constant offset has been added to the top four spec- 
tra for clarity, and the days —5, 0, and +12 spectra have been 
multiplied by a factor of 3 before adding the offset to make the 
emission line visible on this display scale. The emission line ap- 
pears to be stronger at late times primarily because the supernova 
continuum on top of which it is superimposed has faded by more 
than an order of magnitude. 

absorption. 27 In light of these calculations , plus the con- 
firma tion from SN 20071e and SN 1999cl (jBlondin et all 
2009) that the variable Na absorption phenomenon is not 
unique to SN 2006X, we now reconsider the nature of the 
absorbing material. 

4.1. Photoionization and Recombination Timescales 

If the changes in the Na D absorption lines are a result 
of recombination, then the spectra tell us directly that 
the recombination timescale is ~ 10 days (e.g., Fig. [3]). 
Since 



1 



^rccomb — 



( \ It ( 



(i) 



where a, the Na I recombin ation coefficien t, is 1.43 x 
10~ 13 cm 3 s" 1 at T = 10 4 K (jBadnellll200l see g4~2]for 
a justification of this temperature range), we find that 
the electron density must be n e w 8 x 10 6 cm" 3 . (For 
T = 2 x 10 3 K, a = 6.18 x 10" 13 cm 3 s -1 , and the 
required electron density will be lower by a correspond- 
ing factor.) We note that such a high electron density 
necessitates that hydrogen must be partially or mostly 
ionized; no other plausible sour ce could provide so m any 
electrons (as mentioned in § 13.71 Mazz ali et al.ll2005al sug- 
gest that electrons provided by hydrogen can also explain 
the origin of high-velocity features, as were present in 
SN 20071e, at early times). Even with full ionization, 
though, the physical density required in the absorbing 
material is quite high. 

In order for a significant fraction of neutral sodium to 
be present, the photoionization rate should be compara- 

27 Note, however, that the underprediction of Na I is a known 
issue in ionization models (e.g., Mazzali et al. 1997). The reason for 
this problem is not understood, but it could provide an explanation 
for the observed EW ratio between Na I and Ca II exceeding the 
expected one by a large factor (see £|4.2| I, 



ble to or slower than the recombination rate, or alterna- 
tively, the recombination time s hould be shorter than the 
photoionization time. Following lMurrav et alj (|2007f ). we 
write 



A-KT 2 {hu) 

«Na I-^UV 



(2) 



where r is the distance between the SN and the absorb- 
ing material, QNa i ~ 10~ 19 cir2 is the photoionization 
cross section (|Verner et al.l 119961 ) . Ltjv is the UV lumi- 
nosity of the supernova, and (hu) corresponds to the pho- 
ton energy required to ionize a Na I atom. We do not 
have UV photometry of SN 20071e, but to the degree 
that UV light curv es of SNe la are fairly homogeneous 
(jBrown et al.l l2~009) , it is reasonable to substitute another 
similar SN. We therefore use the Swift li ght curve of t he 
well-observed normal SN la 2007af from lBrown et al"] to 
estimate the UV luminosity. SN 2007af reached a peak 
Vega magnitude of ~ 16.2 in the Swift UVW2 filter a 
few days before the B-band maximum. The UVW2 fil- 
ter has a central wavelength of 1928 A, a full width at 
half-maximum intensity (FWHM) of 657 A, and a photo - 
metric zeropoint of 17.35 ± 0.03 mag ()Poole et al.l l2008). 
We use the distanc e modulus of SN 2007af, 32.06 mag 
(|Simon et al.1 12007) , to deduce an absolute magnitude 
of Mjjvw2 = —15.86, which suggests an apparent peak 
magnitude for SN 20071e of 16.49. With the zeropoint 
of the UVW2 filter, this mag nitude trans l ates t o a Swift 
count rate of 2.2 counts s _1 . iPoole etldl (|2008h provide 
a conversion between the count rate and the flux density 
that depends only slightly on the source spectrum (r 



6.1 x 10 16 ergs cm 



A 1 ), so integrating over the 



filter bandwidth yields a flux of 8.8 x 10 -13 ergs cm" 2 s. 
We therefore estimate a UVW2 luminosity for SN 20071e 
of 9 x 10 40 ergs s" 1 . Note that this bandpass corresponds 
reasonably closely to the wavelengths that can ionize Na I 
atoms most efficiently (A < 2412 A). Using Equation^ 
we can now calculate that the ionization timescale is 



= 1.1 



10 



16 



(3) 



For r < 3 pc, then, Ti on <C r rccom b, and this estimate 
suggests that the Na should be essentially fully ionized 
around maximum light. 

Although the UV luminosity of the SN declines with 
time, the Swift data s how that the declin e is rather slow: 
~1 mag in 10 days (H rown et al J 120091) . Even several 
weeks after peak, when the high-resolution spectroscopy 
demonstrates that significant changes in the Na D line 
profile have already occurred, Equations [T] and [3] indicate 
that recombination (at the assumed densities) should 
only be occurring at relatively large radii, calling into 
question the hypothesis that circumstellar material is 
responsible for the varying Na absorption. These simple 
ca lculations are i n qualitative agreement with the results 
of IChu gai (20Ql) and indicate that more sophisticated 



modeling is required. 

4.2. Photoionization Modeling 

Because the simple ionization and recombination 
timescale arguments do not provide obvious answers as 
to why the Na I column density is changing, we now 
consider more detailed photoionization models of the 
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SN environment. Starting with a synthetic UV spec - 
trum for a Type la supernova from iNugent et al.1 (|1995l ) 
as the ionizing source, we carry o ut calculations with 
the photoionization code Cloudy (|Ferland et al.l 119981 
ver. 07.02). We use the UV spectra of SN 2001eh 
and SN 200 lep (jSauer et al.ll2008l ) to normalize the total 
number of UV-IR photons (1500-12000 A) in the model 
spectrum. 28 The model spectrum is only available for 
the time of maximum brightness, so we assume that the 
light curve at all UV wa velengths follows th e SN 2007af 
UVW2 light curve from iBrown et alj (|2009h . (We have 
to combine observations of multiple SNe here because 
there are few nearby type la supernovae with both good 
UV light curves and UV spectroscopy.) This assumption 
means that the calculations are not explicitly time de- 
pendent, but we do not expect non-equilibrium effects 
to produce qualitative changes in the results discussed 
below. We also note that, while at longer UV wave- 
lengths the model can be compared to observed spectra 
for validation, at the critical far-UV wavelengths that 
ionize hy drogen no such comparison with real d ata is 
possible. Bo rkowski. Blondin. fc Reynolds! (|2009f ) have 
recently carried out very similar calculations to inves- 
tigate the possibility of detecting CSM interactions in 
spectra obtained before maximum light. 

Guided by the hypothesis that recombination is re- 
sponsible for the varying Na I column density, a few iter- 
ations of the model demonstrate that very large hydrogen 
densities are needed to produce significant Na recombina- 
tion (for a solar abundance ratio). We therefore assume 
a hydrogen density of n — 2.4 x 10 7 cm' 3 (somewhat 
higher than the electron density inferred in § 14. 1[) . We 
then calculate the Na I, Ca II, and H II fractions as a 
function of time and distance away from the SN. The re- 
sults are illustrated in Figure |T0j For distances greater 
than ~ 0.1 pc the fraction of Ca II ions is essentially 
constant with time, while the Na I fraction increases by 
approximately an order of magnitude over a time span of 
30 days. Closer to the SN, the Ca becomes more highly 
ionized (and the Ca II fraction varies with time) and the 
amount of neutral Na is negligible, supporting the idea 
that the source of the absorption must be located at dis- 
tances larger than 10 16 cm. 

Over the range of parameter space where these cal- 
culations indicate that the Na I column density will 
vary and the Ca II column density will not, the EW 
of the Na D 2 line is predicted to be less than or equal 
to the EW of the Ca K line. At distances of less than 
1 pc, the Ca absorption should be > 100 times as strong 
as the Na absorption. Yet, we observe that EWn s d 
is larger than EWc a h&k at all epochs (see Tables [2] 
and Q] and Fig. [7]). If the model described above is 
correct, we therefore require that nearly all of the Ca 
atoms in the CSM are locked up in dust gr ains, as is also 
typical in dense interstellar clouds (e.g., ISpitzerl I1954 
Howard. Wentzel. fc McGedll963t fHerbid 119681 : iMortonl 
19751) . In the immediate vicinity of a supernova, the de- 
struction of dust grains is obviously a possibility. How- 

28 The normalization is chosen such that 

47rdg N f Xjf^dX = N SN f \J™ odel d\, where J™ dcl is the 
model spectrum, and JJ N and dg^ are the SN spectrum and 
distance. The normalization factors for SN 2001ep and SN 2001eh 
agree to within 30%, so we use TVgN = 0.5(N ep + N e ^). 



ever, given sublimation temper atures of ~ 1500 K for si l- 
icate and carbonaceous grains (|Draine fc Salpeteri[l979h . 
calculations of grain temperatures w i th Clo udy using the 
dust grain model of Ivan Hoof et all (2004) indicate that 
dust can survive at distances greater than 0.05 pc (un- 
shaded region in Fig. I10p over the relevant time period. 

4.3. Physical Conditions in the Absorbing Gas 

The results presented in § 14.21 point to distances of 
~ 0.1 pc for the absorbing material. Figure [10] demon- 
strates that the EW of the Ca II absorption would re- 
main constant for distances as small as ~ 0.02 pc, but 
the fraction of neutral sodium so close to the supernova 
is negligible. Distances greater than 1 pc are also al- 
lowed by the photoionization modeling, but then the 
ionization fraction is very low (as a result of the weak 
UV field), requiring extremely high physical densities to 
produce large numbers of electrons. Of course, at such 
distances the galactic UV field may contribute signifi- 
cantly to the photoionization rate, rendering our cal- 
culations no longer applicable. At d = 0.1 pc, the 
gas temperature is ~ 5000 K, and densities of «h ~ 
2.4 x 10 7 cm -3 are required to produce significant Na 
recombination. With such a high density, substantial 
amounts of dust could be present, and strong depletion 
of Ca atoms onto dust grains would explain the large 
observed ratio of Na I to Ca II. Given the derived den- 
sity and a solar Na abundance of 12 + log(Na/H)= 6.17 
(|Asplund. Grevesse. fc Sauva]||2005[ ). the total density of 
Na atoms is riNa = 35 cm -3 . The Na I fraction un- 
der these conditions is ~ 5 x 10 -4 , yielding n-^ a i = 
1.8 x 10~ 2 cm -3 . In order to obtain a column density of 
~ 2.5 x 10 12 cm -2 , as derived in § 13.41 the path length 
through the absorbing cloud must be ~ 1.4 x 10 14 cm. If 
such a clump of material were roughly spherical its mass 
could be quite small (~ 2 x 10~ 7 M Q ), but a clump of ra- 
dius 10 14 cm at a distance of 10 17 cm would only occult a 
small fraction of the SN photosphere (radius ~ 10 15 cm) 
after maximum light and thus could not be responsible 
for an absorption line as deep as that observed in Fig- 
ure [3] In order to produce a covering fraction closer 
to unity, there must either be large numbers of small 
clumps, or perhaps a larger sheet of material with per- 
pendicular dimensions of ~ 10 cm (the rather high den- 
sities involved prob ably rule out the f ull spherical shell 
model proposed by iPatat et~ai]|2007al ). The total CSM 
mass along the line of sight (presumably more CSM is 
present in other directions) would then be in the range 
of 10~ 5 — 10~ 2 M Q depending on the geometry. 

Is it plausible for such a cloud (or clouds) to exist in 
the SN progenitor system? The measured velocity of 
the variable absorbing component is 2123 km s _1 , while 
the local host galaxy ISM velocity near the position of 
the SN is rs 2130 km s _1 . Assuming that the velocity 
measured for the local ISM is representative of the pro- 
genitor system, then the variable absorption component 
has a line-of-sight velocity of ~ 10 km s _1 relative to the 
progenitor itself. At a constant velocity of 10 km s" 1 , 
material starting near the center of the progenitor sys- 
tem would take ~ 10 4 yr to reach a radius of 0.1 pc, 
but factoring in the gravitational deceleration the actual 
travel time would be ~ 3000 yr if the material originated 
in the red giant wind or ~ 10 yr if it was blown off from 
the white dwarf (which requires a much higher initial ve- 
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Fig. 10. — Results of photoionization calculations performed with Cloudy for a total density of n = 2.4 X 10 7 cm" 3 12 days after 
maximum brightness, (a) Ion fractions for H, Na, and Ca as a function of ionization parameter (or distance along the top axis). The solid 
line shows the Na I fraction, the dashed line shows the Ca II fraction, and the dot-dashed line shows the H II fraction. The black arrow 
illustrates the magnitude of the change in ionization parameter from maximum light to 30 days later, and the shaded gray region on the 
right side of the plot represents the range of ionization parameters for which the dust grain equilibrium temperature is larger than the dust 
sublimation temperature, so that dust grains should not survive (and thus Ca atoms should all be in the gas phase). (6) Temperature and 
equivalent width ratio as a function of ionization parameter. The dotted line represents the ratio of the equivalent width of the Na I D2 
line to that of the Ca II K line (axis scale on the right), and the solid line shows the temperature (axis scale on the left). The Cloudy model 
demonstrates that variable Na absorption, constant Ca absorption, and highly depleted Ca abundances can be obtained for distances larger 
than ~ 0.1 pc away from the site of the explosion. 



locity). Regardless, this timescale is certainly far shorter 
than the age of the progenitor system or its evolutionary 
timescale, so producing CSM at such distances does not 
appear to be problematic. The survival of dense clumps 
similar to what we are proposing here, though, may be an 
issue: the thermal pressure in such a cloud is quite large, 
suggesting that the cloud should expand and dissipate 
on timescales of decades unless the average CSM pres- 
sure is also very high. We do not currently have enough 
information to determine whether such clumps of mate- 
rial could form from the wind of the companion alone, 
the interaction of the wind with surrounding material, or 
a nova eruption on the surface of the white dwarf. 

5. SUMMARY AND CONCLUSIONS 

We have presented a series of high-resolution spectra 
of the nearby Type la supernova SN 20071e, a normal 
SN la in the subgroup that exhibits broad lines and high 
velocity gradients, similar to SN 2002bo and SN 2006X. 
The spectra reveal a complex set of absorption systems in 
the Na D and Ca H & K lines. Our observations demon- 
strate that one of the Na absorption components changes 
systematically with time, increasing its equivalent width 
by 85 mA from maximum light to three months later. 
SN 20071e is the third SN la in which this effect has been 
seen, after SN 2006X (|Patat et al.ll2007al) and SN 1999cl 
(jBlondin et al.ll2009f) . As in SN 2006X, we detect no vari- 
ations in the Ca II absorption column over the course of 
our observations. 

In the initial study of SN 2006X, iPatat et all (|2007al) 
modeled the varying absorption as occurring in several 
shells of circumstellar gas resulting either from the stel- 
lar wind of the companion to the supernova progenitor 



(presumed to be a re d giant) or rec urrent nova erup- 
tions. Subsequently, iChugail (|2008l ) argued that the 
red giant wind scenario could only reproduce the ob- 
serv ed behavior wi th an unrealistically high wind den- 
sity. Chugai (2008) did note that a red giant wind plus 
a dense clump or shell of material might be compatible 
with the observed Na I and Ca II optical depths, but 
concluded that this solution was unlikely because of the 
difficulty of creating a sufficiently dense shell. Moreover, 
both SN 2006X and SN 1999cl were extremely heavily 
reddened, suggesting the possibility of a connection be- 
tween variable absorption and properties of the inter- 
stellar medium along the line of sight. The much lower 
extinction toward SN 20071e demonstrates that variable 
absorption also occurs in supernovae that are not seen 
through an unusually dense, dusty ISM. Conversely, it is 
now also cle ar that some h i ghly obscured objects such as 
SN 2003cg dBlondin et alJl2009h and SN 2008fp (Patat 
et al., in preparation) do not exhibit variable absorption. 
We therefore carried out further photoionization calcula- 
tions to investigate whether circumstellar material could 
explain the varying Na and constant Ca absorption pro- 
files. Using a synthetic UV spectrum for a SN la, we 
found that a small (10 14 cm), dense (10 7 cm" 3 ) clump 
or clumps of material located ~ 0.1 pc away from the 
explosion can provide the necessary column density and 
ionization changes with time to account for the observa- 
tions. In order to reproduce the observed Na I to Ca II 
EW ratio, this model requires a very large fraction of 
the Ca atoms to be depleted onto dust grains. The most 
likely origin for gas so close to the supernova is that it 
was produced in the progenitor system. If this explana- 
tion is correct, the results indicate a single-degenerate 
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progenitor system for SN 20071e, in accord with previous 
suggestions for SN 2006X. 

Through the work of lBlondin et alj (|2009D and our own 
high-resolution spectroscopy of additional SNe la (to be 
described in future papers), it is clear that the fraction 
of SN la events that show the variable Na signatu re is 
relatively small (as predicted bv lPatat et al.l l2007a.b for 
the recurrent nova scenario) . In light of this rarity, it is 
noteworthy that all three of the SNe displaying variable 
absorption share common features in their early-time and 
maximum-light spectr a, classifying them in the broad 
line group according to lBranch et al.1 (120091) an d the high 
velocity gradient group of iBenetti et al.l (|2005l although 
the HVG classification of SN 20071e is less firm). With 
these groups comprising ~ 30% of SNe la, the odds of 
three SNe exhibiting both variable absorption and broad 
lines/high velocity gradients by chance (i.e., if there is 
no physical connection between the two) is low. Our 
observations of SN 20071e therefore support earlier spec- 
ulations of a relationship between high velocity gradients 
and variable Na absorption. 

These results still leave open the three explanations for 
variable absorption proposed in the literature: (1) SNe 
with high velocity gradients and variable absorption oc- 
cur in a distinct progenitor system or with a different 
explosion mechanism than other SNe la; (2) SNe with 
high velocity gradients and variable absorption occur in 
the same progenitor systems as other SNe la, but with 
a peculiar geometry that yields higher observed veloc- 
ities and makes variable absorption more likely; or (3) 
variable absorption is simply a line-of-sight effect that 
is independent of the properties of the supernova itself. 
If the variable absorption originates in the ISM, then it 
is much more likely to be observed for objects that are 
seen through large ISM columns. The low extinction of 
SN 20071e, suggesting that it is not obscured by a sub- 
stantial amount of ISM material, renders the third possi- 
bility less likely. Further high-resolution spectroscopy of 
BL/HVG SNe will be needed to continue investigating 
the first two scenarios and reveal whether high veloc- 
ity gradients are a sufficient condition for variable Na 
absorption or merely a necessary one. It is also worth 
noting that in all three SNe with variable absorption, 
the absorption lines are systematically blueshifted with 
respect to the local ISM velocity, further pointing to a 
connection between the absorbing material and the su- 
pernovae themselves. 

In addition to variable Na absorption, we also de- 
tect Ha emission in the spectrum of SN 20071e, possi- 
bly originating from the supernova itself. The narrow 
Ha emission line (FWHM = 25 km s _1 ) is spatially un- 
resolved, coincident with the position of the supernova, 
and at a slightly different velocity than the surround- 
ing diffuse interstellar emission. The amount of hydro- 
gen that would need to be present in the circumstcllar 
environment to produce the observed emission luminos- 
ity is ~ 10~ 2 — 10 -3 Mq. This feature is present at 
all high-resolution spectroscopic epochs, but as a result 
of its small equivalent width is not visible in any of the 
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low-resolution spectra. The emission-line flux is constant 
within the uncertainties through 12 days after maximum 
light, but possibly strengthens in the observations ob- 
tained several months later. In the best of our spectra, 
H/3, [N II] , and [S II] lines are also detected, but interest- 
ingly [O III] is not, suggesting a very low excitation for 
the ionized gas. If the increase in the Ha flux with time 
is real, then the emitting gas must be associated with 
the supernova, but further observations at late times af- 
ter the SN has faded will be needed to confirm or reject 
this possibility. 
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TABLE 5 
Emission-Line Widths and Fluxes 



03 



Epoch 


Ha EW 


Ha velocity 


Ha line width 


Ha flux 




H/3 flux 


[O III] flux 


[N II] flux" 


[S II] flux b 




[mA] 


[km s- 1 ] 


[km s- 1 ] 


[ergs cm~ 2 s~ 1 


] 


[ergs cm -2 s -1 ] 


[ergs cm~ 2 s _1 ] 


[ergs cm~ 2 s _1 ] 


[ergs cm -2 s _1 ] 


day —5 


22 ±3 


2137.8 ± 1.1 


25.9 + 4.5 


(1.5 + 0.2) X 10" 


16 


c 


< 5.8 x 10~ 17 






day 


14 ±4 


2140.1 ± 1.7 


17.9 + 4.0 


(1.2 + 0.3) X 10" 


16 










day +12 


28 ±5 


2136.4 ± 1.4 


25.4 + 3.3 


(1.5 + 0.3) x 10" 


16 










day +84 


288 ± 9 


2139.3 ± 0.3 


30.8 + 0.7 


(2.1 + 0.1) x 10" 


16 


3.9 + 0.7 x 10~ 17 


< 1.5 x 10~ 17 


(6.0 + 0.7) x 10~ 17 


(2.7 + 0.5) X 10~ 17 


day +90 


301 ± 41 


2139.5 ± 1.0 


24.5 + 2.6 


(1.8 + 0.3) x 10" 


16 







a The listed value is for the stronger [N II] A6583 line. The listed value is for the stronger [S II] A6717 linc. c We have listed upper limits only where the data provide useful constraints; [O III] 
upper limits are at 4 a significance 
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